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Abstract 

This Letter presents a search for tb resonances in 1.04 fb 1 of LHC proton-proton collision data 
collected by the ATLAS detector at a center-of-mass energy of 7 TeV. Events with a lepton, missing 
transverse momentum, and two jets are selected and the invariant mass of the corresponding final 
state is reconstructed. The search exploits the shape of the tb invariant mass distribution compared 
to the expected Standard Model backgrounds. The model of a right-handed W' R with Standard Model- 
like couplings is chosen as the benchmark model for this search. No statistically significant excess 
of events is observed in data, and upper limits on the cross section times the branching ratio of 
W' R resonances at 95% CL lie in the range 6.1-1 .0 pb for W' R masses ranging from 0.5 to 2.0 TeV. 
These limits are translated into a lower bound on the allowed right-handed W' R mass, giving m W ' R > 
1.13 TeV at 95% CL. 
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This Letter presents a search for tb resonances in 1.04 fb _1 of LHC proton-proton collision data 
collected by the ATLAS detector at a center-of-mass energy of 7 TeV. Events with a lepton, missing 
transverse momentum, and two jets are selected and the invariant mass of the corresponding final 
state is reconstructed. The search exploits the shape of the tb invariant mass distribution compared 
to the expected Standard Model backgrounds. The model of a right-handed W' R with Standard 
Model-like couplings is chosen as the benchmark model for this search. No statistically significant 
excess of events is observed in data, and upper limits on the cross section times the branching ratio 
of W' R resonances at 95% CL lie in the range 6.1-1.0 pb for W' R masses ranging from 0.5 to 2.0 TeV. 
These limits are translated into a lower bound on the allowed right-handed W' R mass, giving m w i > 
1.13 TeV at 95% CL. 
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This Letter presents a search for tb (tb or tb) resonances 
using data collected in 2011 by the ATLAS detector [1] 
at the Large Hadron Collider (LHC), corresponding to 
an integrated luminosity of 1.04 ± 0.04 fb _1 [2j, [3| from 
pp collisions at a center-of-mass energy of 7 TeV. These 
resonances include new heavy gauge bosons such as the 
W' boson. The W boson is a charged heavy gauge bo- 
son that is predicted in many extensions of the Standard 
Model (SM) such as universal extra dimensions Q and 
little Higgs models [5(. If the W boson is assumed to 
have similar coupling strengths to those of the SM W 
boson, searches in the W' — > tv decay channel, where 
£ is a charged lepton, are the most sensitive. However, 
the W —> tb channel is competitive if W — > Iv decay 
is suppressed. For example, for a right-handed W' R this 
can happen if the right-handed neutrino, v R , is heavy 
enough to prevent W' R — > lv R decay [fj. The model of 
a right-handed W' R with SM-like couplings is chosen as 
the benchmark model for the analysis presented in this 
Letter. The W' R — > tb decay channel has been searched 
for at the Tevatron 0, The best previous limit on 
a W' R with Standard Model-like couplings of the W to 
quarks was set by the DO experiment and excludes a W' R 
mass below 890 GeV at 95% confidence level. 

The innermost part of the ATLAS detector [§] , a track- 
ing system in a 2 T axial magnetic field, measures the 
momentum of the charged particles produced in the colli- 
sions. Outside of the solenoid are the calorimeter subsys- 
tems, which measure the electron, photon, and hadronic 
particle energies, and the muon spectrometer which is 
used to identify and measure the momentum of muons 
in a toroidal magnetic field. A three-level trigger sys- 
tem [l(| reduces the event rate and selects the events for 
analysis. 

The tb resonances are searched for in the tb — » Ivbb 
decay channel, where the lepton, £, is either an electron 
or a muon. W' R signal events are simulated to leading 
order (LO) with the PYTHIA v6.421 [TO] Monte Carlo 



(MC) generator, using the MRST2007 LO* parton distri- 
bution functions (PDFs) [13]. Seven signal samples are 
simulated, with different W' R mass assumptions, ranging 
from 500 GeV to 2.0 TeV, as reported in Table |H The 
respective signal cross section times the branching ratio 
values are computed at next-to- leading-order (NLO) [l3j |. 
using CTEQ6.6 PDFs [l4|. 



[GeV] 



B{W' R -> tb) 



ax B [pb] 



500 
750 
1000 
1250 
1500 
1750 
2000 



0.298 ± 0.002 
0.319 ±0.001 
0.326 ± 0.001 
0.328 ± < 0.001 
0.330 ± < 0.001 
0.331 ± < 0.001 
0.332 ± < 0.001 



54.6 ±2.1 
10.9 ±0.6 
2.92 ±0.18 
0.91 ±0.07 
0.31 ±0.03 
0.11 ±0.01 
0.044 ±0.005 



TABLE I: NLO branching ratios, B{W' R -> tb), and W' R 
production cross section times the branching ratio value, 
a(pp — s> W' R ) x B(W' R — !■ tb), in pp collisions at 7 TeV center- 
of-mass energy [L|| • The uncertainties on the branching ratios 
are due to the top quark mass uncertainty. The uncertainties 
on the cross sections include statistical, a s , NLO renormal- 
ization and factorization scales, and PDF uncertainties. 



Data-driven methods and MC simulated samples are 
used to estimate and model backgrounds. The ti pro- 

SET 



cess is simulated with the MC@NLO v3.41 \m,UM MC 
generator, assuming a top quark mass of 172.5 GeV, and 
using the CTEQ6.6 PDFs. The parton shower is added 
using the HERWIG and JIMMY MC genera- 
tors. The tt cross section is obtained from the approx- 
imate NNLO prediction calculated with the HATHOR 
program using the MSTW2008 NNLO PDF sets ^ . 
The single top quark processes are simulated using the 
AcerMC v3.7 [21j MC generator and hadronization is 
performed with the PYTHIA MC generator; the cross 
section is calculated to approximate NNLO [22H24] using 
the CTEQ6.6 PDFs. Diboson processes are simulated 



2 



using the HERWIG v6.5 MC generator and their cross 
sections are obtained at NLO using the MCFM 25[ pro- 
gram with the MSTW2008 PDFs. The MC samples sim- 
ulated with the AcerMC and HERWIG MC generators 
use the MRST2007 LO* PDFs. Vector boson produc- 
tion in association with jets (W+ light jets, Wbb, Wcc, 
Wc and Z+jets with up to five additional partons) is 
simulated using the ALPGEN v2.13 26] MC generator, 
coupled with the CTEQ6L1 PDFs [lj] and hadroniza- 
tion is performed with the HERWIG and JIMMY MC 
generators. In these samples, additional jets can be cre- 
ated from the parton shower. In order to avoid double 
counting between the inclusive W+n parton samples and 
the parton shower, overlaps are removed following the 
MLM matching prescription [26] . A cross section correc- 
tion factor is applied to the LO W / 'Z+jet cross sections 
computed by comparing the LO and NLO predictions 
from the FEWZ [27J program. The Wc cross section 
correction factor is obtained using the MCFM [28| pro- 
gram with the CTEQ6.6 PDFs. All samples are passed 
through the full simulation of the ATLAS detector ^ 
based on GEANT4 [3(J and are then reconstructed us- 
ing the same procedure as collision data. The simulated 
samples include the effect of multiple pp collisions per 
bunch crossing (pile-up) which on average is six events 
per bunch crossing. In order to ensure a good description 
of the energy scale and resolution, the trigger, the recon- 
struction and identification efficiency, corrections based 
on comparisons between data and MC events are applied 
to the simulated signal and background samples. The 
corresponding scale factors are obtained as a function of 
the object kinematics, resulting in final corrections of the 
order of a few percent. 

Candidate events are identified using single high trans- 
verse momentum electron and muon triggers and strin- 
gent detector and data quality requirements. For each 
candidate, two jets, one isolated charged lepton, and 
missing transverse momentum, E^ 1BS , are required. The 
definition of the objects and details of a similar event se- 
lection including lepton isolation requirements are given 
in Ref . 
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The reconstructed charged lepton is required 
to have a transverse momentum p^ > 25 GeV to ensure 
a constant trigger efficiency, |f?| < 2.5 for a muon 32- 
[3] and 1 77 1 < 2.47 for an electron [35| (the calorimeter 
transition region 1.37 < |ry| < 1.52 is excluded), and to 
lie within A_R = y/{Arj)^ + (A0) 2 < 0.15 of the corre- 
sponding triggered lepton. Jets are reconstructed from 
energy clusters in the calorimeters with the anti-fct al- 
gorithm [3(| with a radius parameter R = 0.4 and cal- 
ibrated to the hadronic energy scale 37J. Exactly two 
jets with pt > 25 GeV and |?7| < 2.5 are required in the 
event, and at least one of them must be tagged as a b- 
jct. The 6-tagging algorithm uses measurements of the 
impact parameters of tracks and the properties of recon- 
structed vertices; these are combined in a neural network 
to extract a tagging decision for each jet [HI . Based on a 



ti MC sample, the working point is chosen at a 6-tagging 
efficiency of 57%, leading to a light-quark tagging proba- 
bility of 0.2% derived from the same sample. To account 
for the differences between observed and simulated jet 
Pt and 77 distributions, the 6-tagging efficiency and the 
corresponding scaling factors to be applied to MC are 
derived from data [39(. Events before applying any b- 
tagging are referred to as pre-tagged events. Events where 
one or both jets are 6-tagged are referred to as single- or 
double-tagged events respectively. 

The E™ lss is calculated using calorimeter energy clus- 
ters [13] calibrated according to the reconstructed physics 
object to which they are associated 4l|; events are re- 
quired to satisfy E™ ms > 25 GeV. The background contri- 
bution from multiple hadron jets (multijet background) 
is reduced by imposing a requirement on the sum of 
the W boson transverse mass rriT(W) [42[ and E^ 1BS '- 
m T (W) + E^ iss > 60 GeV 43]. After applying all se- 
lection criteria, the acceptance times efficiency for W' R 
signal events with mw> =1.0 TeV is 1.38% for single- 
tagged events and 0.49% for double-tagged events. 

One of the most important backgrounds for the tb res- 
onance search comes from W production in association 
with either heavy-flavor jets, or light-flavor jets misiden- 
tified as 6-jets. Multijet production is another source of 
background, when either a hadronic jet is misidentified 
as a lepton, or when a real high-px lepton from scmilcp- 
tonic decay of a heavy hadron within a jet fulfills the 
selection requirements. Another important background 
comes from ti pair production in the case that one W 
boson decays leptonically and the decay products of the 
other W boson are lost due to the detector acceptance. 
Other smaller backgrounds come from single top produc- 
tion, diboson production and Z+jet events. 

Kinematic variable distributions for the W"+jet back- 
ground are taken from MC samples, while the overall nor- 
malization and flavor composition are derived from data; 
this is done after rejecting signal-like events with the tb 
invariant mass, m t b which is described later, satisfying 
m t b > 500 GeV. In each jet multiplicity bin, the num- 
ber of Vy+jet events in the data is assumed to be the 
difference between the number of observed data events 
and the number of events estimated for SM non-Vy+jet 
processes including the multijet process estimated from 
a data-driven method. The overall VF+jet normalization 
factor is the ratio of the number of W+jet events in the 
data to the number of VF+jet events in simulation. The 
flavor composition of the VF+jet background is estimated 
by comparing the MC prediction to data while its depen- 
dence on jet and 6-tagging multiplicity is modeled using 
MC. The fractions of Wbb, Wcc, Wc and PF+light jet 
components of the total W+jet MC are scaled such that 
the background sum equals the observed data in three 
separate samples: a single-tagged one-jet sample and the 
pre-tagged and single-tagged two-jet samples. The same 
scale factor is used for Wbb and Wcc. 
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The multijet background normalization and the shape 
of each distribution are obtained from data. The shape 
of each multijet background distribution is taken from a 
data sample which requires a jet instead of an isolated 
lepton. This jet is required to have a detector signature 
similar to an electron: it must have pr > 25 GeV and be- 
tween 80% and 95% of its energy deposited in the electro- 
magnetic section of the calorimeter. The jet must also be 
associated with at least four tracks. The normalization is 
estimated using a binned likelihood fit to the _E™ ISS distri- 
bution in data in which the normalization of the VT±jet 
and the multijet components is allowed to vary. The 
fit is performed separately in the pre-tagged, single- and 
double-tagged samples, after applying all selection crite- 
ria except the E™ 1SS cut. The uncertainty on the multijet 
rate is 50% for pre-tagged and single-tagged events, while 
it amounts to 100% for double-tagged events. The un- 
certainty is estimated by using the m^(W) distribution 
instead of the E™ lss distribution in the binned likelihood 
fit, and by using multijet background models built from 
data samples with low and high numbers of pp collisions 
per event. 

The tt, single top, Z+jet and diboson events are nor- 
malized to the theoretical cross sections and the shape of 
each distribution is taken from the MC simulation. 

Based on the theoretical predictions shown in Table |H 
the numbers of single- and double-tagged W' R signal 
events expected in 1.04 fb _1 are listed in Table UI1 as 
a function of mw' R ■ Table IIIII lists the expected back- 
ground yields. 



m w i [GeV] Single-tagged Double-tagged 
ii 

500 973 ± 37 455 ±17 

750 174 ± 9 77 ±4 

1000 42 ±3 15 ±1 

1250 11 ±1 3.9 ±0.3 

1500 3.2 ±0.3 1.0 ±0.1 

1750 1.0 ±0.1 0.26 ±0.03 

2000 0.36 ± 0.04 0.09 ±0.01 



TABLE II: Predicted signal event yields derived using the 
theoretical cross section times the branching ratio values for 
W' R —¥ tb, for single- and double-tagged two-jet events in 
1.04 fb _1 of data. The uncertainties correspond to the NLO 
cross section uncertainties [l3| ]. 

The tb invariant mass is used as the observable to dis- 
criminate signal from background. The neutrino momen- 
tum in the decay tb — > Ivbb is computed assuming the 
transverse component to be equal to E^ lss , and extract- 
ing the longitudinal component (p z ) by constraining the 
l-v invariant mass to mw = 80.42 GeV. This gives a 
quadratic equation in p z and the solution with the smaller 
\p z | is used. If the solution is complex, only the real part 
is taken and the imaginary part is neglected. 

Figure [1] shows the data and expected background dis- 



Samples 


Single-tagged 


Double-tagged 


W+ jets 


5970 ± 1000 


290 ±180 


Multijets 


1120 ±560 


47 ±47 


tt 


1560 ± 130 


o/?rv 1 on 

360 ± 30 


Single top 


1240 ± 90 


120 ±10 


Diboson, Z+jets 


320 ± 120 


14 ±2 


Total prediction 


10200 ± 1200 


830 ± 190 


Data 


10428 


844 



TABLE III: Predicted background event yields compared to 
the total observed event yields for single- and double-tagged 
two-jet events in 1.04 fb~ of data. All VK±jet samples are 
scaled by the factors determined from data, with the uncer- 
tainties also derived from data. The multijet estimation is 
from the fitting method with a 50% (100%) uncertainty for 
single(double)-tagged events. All the other predictions are 
derived using the theoretical cross sections and uncertainties. 



tributions of mtb for single- and double-tagged two-jet 
events. The data event with the highest mtb value corre- 
sponds to a single-tagged event with m t b — 2.0 TeV. The 
Bumphunter tool [44j is used to search for a local ex- 
cess in the data due to the production of a tb resonance. 
This tool is used to test the consistency of the data with 
the SM background only hypothesis, comparing the data 
to the SM prediction over the spectrum of the tb invari- 
ant mass, scanning over sliding mass windows from 0.5 
to 2.0 TeV. The width of the mass windows is chosen 
to be constant in log(mtb) as shown in Fig. [1] to deal 
with low background MC statistics in the higher mass 
bins. This comparison has been performed for single- 
and double-tagged events separately. The region with the 
highest data-background difference is 1024-1129 (764- 
842) GeV for single (double)-tagged events. The probabil- 
ity of observing the SM background fluctuating up to or 
above the number of observed data events in these regions 
is 0.66 for single-tagged events and 0.72 for double-tagged 
events. These values, which are based on the statistical 
error only, indicate that there is no significant evidence 
for tb resonances in the observed data. 

Systematic uncertainties from various sources affecting 
the background and the signal acceptance (rate uncer- 
tainty), as well as shape changes in the invariant mass 
distribution (shape uncertainty) are considered. 

The jet energy scale and the uncertainty on the b- 
tagging scale factors are the dominant systematic un- 
certainties for the signal. The background normalization 
yields are the dominant systematic uncertainty for the 
background contribution. The jet energy scale uncer- 
tainty is evaluated by scaling la up or down the energy 
of each jet. The 6-tagging scale factors are px-dependent 
and have an uncertainty between 8% and 20%. The 
multijet background uncertainty has already been de- 
scribed. The uncertainty on the normalization of the 
VK+jets background and its flavor composition include 
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FIG. 1: The distribution of mtb for single-tagged (top) and 
double-tagged (bottom) two-jet events in data compared to 
Standard Model expectations. The expected W' R signal, nor- 
malized to the theoretical cross section times the B(Wr — > tb) 
values from Tabled) has been scaled by a factor of 10. The ef- 
fect of the jet energy scale (JES) uncertainty on the predicted 
background is shown, as are the data statistical uncertainties. 
The bin width is constant in log(m t j,). The highest bin in each 
plot includes overflows. 



both systematic contributions and a statistical contribu- 
tion from the limited size of the sample. The W+jet fla- 
vor uncertainties are treated as fully correlated between 
Wbb and Wcc and uncorrelated otherwise. Theoretical 
cross section uncertainties for the top (tt and single top) , 
diboson and Z+jet backgrounds of 10%, 5% and 60% 
are assigned, respectively. The Z+jet theoretical cross 
section uncertainty is estimated based on the variation 
of ALPGEN parameters, and a relative uncertainty of 
50% on the heavy-quark contributions, but it has a very 
small impact on the result due to the small contribu- 
tion of Z+jet events. Systematic uncertainties due to 
the residual differences between data and MC simulation 
for the reconstruction and energy calibration of jets, elec- 



trons and muons are estimated to have a small impact on 
the result. The uncertainty on the integrated luminosity 
is 3.7% [U. The uncertainty on the background model- 
ing in the mtb distribution is evaluated using pre-tagged 
data and found to be negligible. 

An uncertainty due to the MC event generator is esti- 
mated by comparing MCS 
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6NLO and POWHEG 
for tt and AcerMC and MC@NLO for single top events. 
The uncertainty in parton shower modeling is estimated 
by comparing two POWHEG tt samples for which the 
hadronization is performed by PYTHIA or HERWIG. 
Uncertainties from modeling the amount of initial and 
final-state QCD radiation are also taken into account. 
The uncertainty due to the specific choice of PDFs 
in the simulated events is determined by re-weighting 
the MC events using the NNPDF20, MSTW2008 and 
CTEQ6.6 Q eigenvector PDF sets. Finally, an uncer- 
tainty to account for the limited MC sample sizes is also 
included. 

No significant data excess is identified for any value of 
mtb, and an upper limit on the W' R — ► tb production cross 
section (a) times the B(W R —> tb) at 95% credibility- 
level (CL) is determined using a Bayesian approach as- 
suming flat priors 47J. The likelihood function used is 
the product of the Poisson probabilities over all mass 
bins 48| per channel. The combination of single- and 
double-tagged events is done by extending the likelihood 
function; the joint likelihood is the product of Poisson 
probabilities for each individual bin in each channel. Sys- 
tematic and statistical uncertainties are incorporated and 
treated as nuisance parameters with a Gaussian proba- 
bility density function. Figure [5] shows the observed and 
the expected limits from single- and double-tagged events 
combined. Observed (expected) upper limits obtained on 
a(pp -> W' R ) x B{W' R -> tb) at 95% CL lie in the range 
6.1 - 1.0 (4.5 - 1.4) pb for W' R masses ranging from 0.5 
to 2.0 TeV. These a x B limits are also applicable to 
a left-handed W. The a x B limits are converted into 
mass limits using the intersection between the theoreti- 
cal a x B curve as a function of mw' R an d the expected 
and observed a x B limit curves. The corresponding ob- 
served (expected) 95% CL lower limit is: m^i > 1.13 
(1.13) TeV. These are currently the most stringent direct 
limits on production of W' R — > tb. 
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boson mass. The theory curve is also shown. 
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